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ABSTRACT 

Here we present an analysis of the global-integrated mid-infrared emission flux of the Earth based 
on data derived from satellite measurements. We have studied the photometric annual, seasonal, 
and rotational variability of the thermal emission of the Earth to determine which properties can be 
inferred from the point-like signal. We find that the analysis of the time series allows us to determine 
the 24 hr rotational period of the planet for most observing geometries, due to large warm and cold 
areas, identified with geographic features, which appear consecutively in the observer's planetary 
view. However, the effects of global-scale meteorology can effectively mask the rotation for several 
days at a time. We also find that orbital time series exhibit a seasonal modulation, whose amplitude 
depends strongly on the latitude of the observer but weakly on its ecliptic longitude. As no systematic 
difference of brightness temperature is found between the dayside and nightside, the phase variations 
of the Earth in the infrared range are negligible. Finally, we also conclude that the phase variation of 
a spatially unresolved Earth-Moon system is dominated by the lunar signal. 

Subject headings: Earth — infrared: planetary systems — planets and systems: individual(Earth) — 
techniques: photometric 



1. INTRODUCTION 

Finding and characterizing habitable exoplanets is one 
of the key objectives of 21st century astrophysics. Tran- 
sit spectroscopy already allows us to constrain some of 
the at mospheric propertie s of hot exoplanets as compo- 
sition fTinetti et al.|| 2007[ ), tempera ture structure (e.g ., 
Knutson et al.[ ( |2008p joTcirculation ( |Snellen et al.|2010[ ). 
The phase variations of the thermal emission or reflected 
light of hot exoplanets can be extracted from combined 
star+planet photometry, providing another powerful tool 
to characterize the climate of exoplanets. This has been 
achiev ed in transiting configurations (e.g., Knutson et al. 



(|2007[)), in cluding one terrestrial-mass planet (Batalha 
et al.|2' 



et al.|20lTj), but als o for non-transiting hot Jupiters (e.g., 
Cowan et al. (2007)). In theory, the thermal orbital light 
curves of non-transiting short-period rocky planets can 
also be measured with combined lig ht photometry and 



used to detect atmospheric species (Selsis et al. 20111 



or to constrain the radius, albedo and inclination of the 



planet (Maurin et al. 2012). Transit spectroscopy and or- 
bital light curve measurements may be achievable with 
the James Webb Space Te lescope or the Exopl anet Char- 
acterization Observatory (Tinetti et al. 2011) for a few 
habitable planets transiting M stars that can hopefully 



be found in the solar vicinity ( [Belu et al |2011| |Palle 



et al.| 2011 ). The characterization of an Earth analog 
around a G star by transit spectroscopy is much more 
challenging for several reasons: G stars are less numer- 
ous than M stars, the transit probability in the habitable 
zone is 10 times lower than it is for M stars, the orbital 
period (and thus the duration between two transits) is 
significantly longer for G stars, and the planet-to-star 
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contrast ratio less favorable for secondary eclipse spec 
troscopy. Direct detection thus seems to be necessary to 
study Earth analogs around G stars and several instru- 
ment concepts have been proposed (|Traub et al.l|2006 



Danchi fc Lopez|2007| [Trauger k. Traub|2007[|Cash et a 



2009[ ). Depending on the concept, it is either the light 



scattered in the visible or the infrared emission that can 
be detected. Spectroscopy and photometry can then be 
used to derive some of the planet properties. An im- 
portant step toward these ambitious programs is to de- 
termine what level of characterization could be achieved 
when observing the Earth as a distant pale infrared dot. 
In the optical range, broadband photometry can poten- 
tially give us information about the albedo and the cloud 
cover. The uneven distribution of oceans and continents 
enables the measurement of the 24 hr rotation period by 
the autocorrelation of the signal for several viewing in- 
clination a ngles despite the presence of active weather 
(Palle et al.||2004 2008). Visible and near- i nfrared spec- 
troscq p ical studies have be e n made by, e.g. 
" 2011 



2009|), |Cowanet al 



Robinson et al. 
observing rotational and seasonal 



Hearty et al. 
1 2011 ap and" 



|Livengood et al 

variations ot the Earth spectrum and their influence on 
the detectability of the spectral signatures of habitability 
and life. 

In this paper, we provide an integrated mid-infrared 
(4-50 /im) photometric time series model of the Earth, 
with 3 hr time resolution, constructed from over 22 years 
of available satellite data. From this geographically re- 
solved data set, we derived the disk-integrated photomet- 
ric signal of the Earth, seen as a point-source planet. The 
paper is organized as follows: in Section 2, we describe 
our model, the input data, and the assumptions that are 
considered. In Section 3, we describe the analysis of the 
time series and we discuss the retrieving of the rotation 
period, the variability of the signal, and the case of the 
Earth-Moon system. In Section 4, we summarize our 
conclusions. 
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Fig. 1. — Maps of Earth's outgoing mid-infrared radiation. Average over the period 18:00-21:00 UT of 1987 July 1 (a) and 2001 July 1 
(b) together with the subsolar point and the terminator at the mean time. At that time the brighter regions of the Earth correspond to 
the deserts of Sahara, Arabian Peninsula, Atacama, and Arizona. Average over the months 1987 July (c) and 2001 July (d). Average over 
the years 1987 (e) and 2001 (f). 



2. METHODS 

The Sun outputs most of its energy in the visible 
and near-infrared range of the electromagnetic spectrum 
which is absorbed by the Earth and re-emitted to space 
in the mid-infrared range. The equilibrium temperature 
of the planet can be generally defined as: 
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where F$ is the solar incident flux, A is the terrestrial 
albedo, a is the Stefan-Boltzmann constant, Lg is the 
solar luminosity and a is the Earth-Sun distance. For 
the Earth, this temperature is approximately T eq ~255 K 
and the emission reaches its maximum at 10 fim although 
its spectral distribution is strongly affected by the broad 
absorption bands of atmospheric greenhouse gases. The 
emission that a remote observer would detect, however, 
will be highly dependent on the season, the planetary 
region that is viewed and the cloud distribution at that 



time. Here, we derive the variability and periodicities of 
this integrated flux as it would be seen from a remote 
perspective, using satellite data and a geometrical model 
to build the point-like signal received by the observer. 
The interest of this work lies in the facts that land masses 
can emit more infrared radiation than ocean areas, these 
latter have a bigger thermal inertia that make them more 
resistant to temperature changes and the atmosphere of 
the Earth is optically thin enough in some windows of 
the thermal infrared to observe the influence of surface 
features on the emission at the top of the atmosphere 
(TOA). 



2.1. Data 

To construct the time series of the emitted infrared 
flux, we have used TOA all-sky upward longwave (LW) 
flux integrated over the 4-50 /im wavelength interval. 
These data were obtained from the NASA Langley Re- 
search Center Atmospheric Sciences Data Center, they 
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are part of the NASA/GEWE X SRB Projecf] ( jSuttles 
et al. 111989 Gupta et al. 1992 ) and they cover a 22 year 



period from 1983 to 2005 The GEWEX LW algorithm 
used to create these data (Fu et al. 1997) uses a thermal 



infrared radiative transfer code with input parameters 
derived from the Inte rnational Satellite C loud Climatol- 
ogy Projecf1(ISCCP; |Rossow et al.| ( |l99 6|), the Goddard 
EOS Data Assimilation System level-#] (GEOS-4), the 
Total Ozone Mapping Spectrometer^] (TOMS) archive, 
and the TIROS Operational Vertical Soundef] (TOVS) 
data set. The data have a time resolution of 3 fir over the 
whole globe, and a spatial resolution of l°xl° square cells 
in latitude and longitude. Typical maps of the outgoing 
mid-infrared radiation of the Earth, directly represented 
from the GEWEX data, are shown in Figures llTa) and 
(b). The maps represent the average flux over the pe- 
riod 18:00-21:00 UT for 1987 July 1 and 2001 July 1, 
respectively. Some desert regions are noticeable, such as 
Sahara, Arizona or Atacama. Some cold and humid re- 
gions are remarkable too, such as Indonesian clouds, or 
Antarctica. 

2.2. Geometrical Model 

In order to simulate the infrared observations of the 
Earth seen from any given direction, we have built a ge- 
ometrical model of the emitted radiation of the Earth. 
This model calculates at any time the fraction of the 
planetary disk exposed in the direction of the remote ob- 
server and computes the integrated flux F t, s F Q f, s received 
by the observer at a distance d: 



F ohs (t) = jpYl 
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where is the area of each grid cell, cosflj is the co- 
sine of the angle between the surface normal and the 
observer's line of sight, F{ is the TOA flux given by 
GEWEX ( with a resolution of 3 hr) and Fi/n is the 
specific intensity in the Lambertian approximation. We 
obtain a point-like source emission every 3 hr. 

With the purpose of illustrating several possible planet 
viewing geometries, we have defined our observer posi- 
tion by equatorial coordinates at a radius distance of 
10 pc. The equivalent of the declination angle on the 
planet surface gives us the latitude of the sub-observer 
point, which remains constant along the orbit. The 
position of the planet in the orbit and the position of 
the sub-observer's point with time are derived from the 
ephemerides of the substellar point. The e phemeris data 
are ta ken from the JPL Horizons SysterrP ( Giorgini et al 
1996| . 

In order to understand and interpret the variation of 
the Earth emission in relation with geographic features 
and climates, it is relevant to locate the sub-observer 
point by its geographic coordinates. Indeed, its latitude 
remains constant with time because precession and nu- 
tation are negligible during an orbit, and its longitude 

4 http://www.gewex.org/ 

5 http://isccp.giss.nasa.gov/ 

6 http://daac.gsfc.nasa.gov/ 

7 http://jwocky.gsfc.nasa.gov/ 

8 http:/ /www. ozonelayer.noaa.gov/action/tovs. htm 

9 http:/ /ssd.jpl.nasa.gov/?horizons 
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Fig. 2. — Diagram of the model. We take into account the four 
observers per declination angle at a distance of 10 pc in the planes 
defined by the directions of opposition (O), western quadrature 
(WQ), conjunction(C), and eastern quadrature (EQ) and the rota- 
tional axis. The intersection of the direction vectors (red) with the 
sphere of 10 pc (blue) define their position. The sub-observer point 
is defined by the intersection with the planetary surface (black). 



simply varies as Lit, uj being the rotation rate of the 
Earth (with the exceptions of the polar cases where the 
planetary geometry does not change with time). There- 
fore, we can define the sub-observer point by its latitude 
and by its longitude at a given time. This time is set 
to January 1, 0:00 UT of the year considered. We con- 
sider four different initial longitudes for the sub-observer 
point: the meridian of the subsolar point, the morning 
terminator, the meridian of the antisolar point and the 
evening terminator. We call these initial observing ge- 
ometries conjunction (C), western quadrature (WQ), op- 
position (O) and eastern quadrature (EQ), respectively, 
although these terms should normally apply only to an 
observer located in the ecliptic plane (Figure [2]) . 

It is important to note two facts: the first is that each 
observer sees the planet at a certain local hour during a 
whole rotation period and the second, that observers that 
are placed at the same latitude but different longitudes 
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see the same region of the planet at different local hours 
(if the (O) observer sees a given region of the planet dur- 
ing the winter midnight and the summer noon, the (C) 
observer sees the same region during the winter noon and 
the summer midnight). Furthermore, an observer over a 
polar latitude sees the same hemisphere of the planet 
along the time of observation. In the need of a reference 
for time, we define its "local hour" as the Universal Time 
(UT). In this case, the seasonal variability and the diur- 
nal variability (daily change in temperature in a certain 
region) produce the variation of the signal and not the 
changing of the planetary view (rotational variability), 
as we see in Section 3.2. 

2.3. Limb Darkening 

In order to convert the TOA flux into the disk- 
integrated flux received by a distant observer, we assume 
an isotropic (Lambertian) distribution of specific inten- 
sities at the TOA. In reality, specific intensities are not 
isotropic, producing limb-darkening or -brightening de- 
pending on the thermal structure and composition of the 
atmosphere, as well as the wavelength. This approxima- 
tion gives an exact result only if (1) the atmosphere is 
uniform (its structure and composition is the same at 
all latitudes and longitudes) and (2) the TOA flux has 
been obtained by integrating over all angles the actual 
distribution of specific intensities. If these conditions are 
fulfilled, then the angular integration over the planetary 
disk is equivalent to the angular integration of the out- 
going specific intensities done to compute the local flux 
at the TOA. Although condition (1) is obviously not ful- 
filled, the geographic variations of the atmospheric struc- 
ture and the composition of Earth's atmosphere are not 
steep and, as a consequence, significant departures are 
expected only for peculiar regions (for instance, the Sa- 
hara) when observed at the edge of the planetary disk. 
However, in the case where these regions do not remain 
all the time too far (>70°) from the sub-observer point, 
this effect should not affect the amplitude nor the period- 
icity of the photometric variations but only their shape. 
In order to evaluate the error caused (for a uniform atmo- 
sphere) when the outgoing flux is inferred from one single 
direction of propagation (two-streams approximation), 
we have used the PhoenisFjcode (a one-dimensional line- 
by-line code with sphericalgeometry and full angular in- 
tegration of specific intensities with a 0.5° precision), for 
composition and thermal profiles typical of the Earth. 
The version of Phoenix we have used was specifica lly de- 
veloped for Earth- like atmospheres (Paillct 2006). We 
find the maximum error on the disk-integrated spectrum 
to be about 5% for a monochromatic flux. At most wave- 
lengths, and on the wavelength-integrated flux, the error 
is smaller than a few percents. 

3. TIME SERIES ANALYSIS 

Two examples of the annual time series of Earth's out- 
going mid-infrared radiation are plotted in Figure [3j The 
high frequency variability corresponds to the emitted 
mid-infrared flux due to Earth's rotation super-imposed 
to the seasonal variation during the year. It is readily 
observable from the figure that the amplitude of the ro- 

10 http://www.hs.uni-hamburg.de/EN/For/ThA/phoenix/index. 



TABLE 1 

Photometric Variability of the Integrated Earth 
Mid-infrared Flux (Percentage Over the Mean Value) 



Photometric Variability (%) 



Viewing Angle 


Orbital Amplitude 


Rotational Amplitude 


1987 


2001 


Jan 


Apr 


Jul 


Oct 


90°N (N. Pole) 


20 


19 


1 




2 




45°N (Mid-lat) 


14 


14 


4(6) 


4(4) 


8(5) 


5(5) 


0° (Equator) 


6 


5 


5(8) 


5(4) 


8(5) 


4(5) 


45° S (Mid-lat) 


7 


7 


3(5) 


3(2) 


3(2) 


2(3) 


90°S (S. Pole) 


11 


10 


1 


1 


1 


1 



Note. — Mean amplitude values of the orbital (seasonal) vari- 
ability for the years of 1987 and 2001 and of the rotational vari- 
ability over the months of 2001 January, April, July and October 
for (0) and (C) (in parenthesis) observers situated at different 
latitudes. Each value is calculated as the percentage of the av- 
erage variation over the mean value. Data analysis for different 
years/months retrieves similar results. 



tational variability is larger for an observer in the equa- 
torial plane (black) and decreases toward more poleward 
views (violet, magenta). On the contrary, because the 
obliquity of the Earth is about 23.44° and then the vari- 
ation of the annual insolation is larger at higher latitudes, 
the amplitude of the seasonal variability increases for the 
polar geometries and decreases towards the equatorial 
view. As it is expected, the seasonal variation of the 
northern latitude time series is opposite to the southern 
ones due to the seasonal cycle of solar insolation. For 
observers over the Northern Hemisphere, the time se- 
ries reach a maximum during the boreal summer near 
the beginning of August (day 211). The equatorial view 
follows the same pattern, showing a "northern-biased" 
behavior of the planet. This effect is due to the exis- 
tence of large landmasses in the Northern Hemisphere 
that emit more infrared radiation than the oceans and 
make the northern summer hotter (Section 3.2). Table m 
shows the amplitude values of the seasonal variability of 
the globally integrated mid-infrared flux for the years of 
1987 and 2001 and for five viewing geometries: 0° (Equa- 
tor), 45°N, 45°S, 90°N (North Pole) and 90°S (South 
Pole). The amplitude is given in percentage change over 
the mean annual value. In order to conduct a further 
seasonal study of the mid-infrared Earth emission, four 
representative months (January, April, July and Octo- 
ber) were selected per year. For the sake of simplicity, 
we present in this paper the results of the seasonal study 
of 2001 as the rest of the data set give similar results. 
The values given in the table are the values for the ob- 
servers placed in the opposition plane and in the con- 
junction plane (in parenthesis), as they have the same 
view of the planet at different local hours. We can see 
that the orbital amplitude variation of a northern lati- 
tude is twice the value of the southern equivalent lati- 
tude. The seasonal variability dominates the signal in 
most latitudes except for the equator, where seasonal 
and rotational variability are similar. It is notable that 
whereas the mean temperature varies with the season, 
the rotational variability depends strongly on the local 
hour (Section 3.2). Mid-latitude and equatorial view- 
ing geometries show a pronounced and synchronized ro- 
tational variability (up to an 8% change), whereas the 
polar light curves present almost no variations. 
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Fig. 3. — Time series of the mid-infrared emission flux for the two years of 1987 (top) and 2001 (middle and bottom). The sub-observer's 
point is represented by the latitudes 90°N, 60°N, 45°N, 30°N, 0°, 30°S, 45°S, 60°S, 90°S, and 0° longitude at the initial time (January 1, 
0:00 UT) and the direction planes of opposition (O) (top and middle) and conjunction (C) (bottom). 
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Fig. 4. — Autocorrelation functions of the mid-infrared emission flux from Earth. 2001 January, April, July, and October for the latitudes 
90°N, 60°N, 45°N, 30°N, 0°, 30°S, 45°S, 60°S and 90°S. 
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Fig. 5. — Rotational light curves of the mid-infrared radiation emitted from the Earth. 2001 January (top row) and July (bottom row) 
for the latitudes 90°N, 60°N, 30°N, 0°, 30°S, 60°S, and 90°S latitudes, the colors correspond to local hours, hr (black), 6 hr (blue), 12 
hr (green) and 18 hr (red) and the direction planes of opposition (O), conjunction (C), western quadrature (WQ) and eastern quadrature 
(EQ). 



3.1. Periodicities 

The analysis of the cross-correlation of the time series 
with itself or autocorrelation, shown in Figure HJ is a 
technique that allows to determine the rotation period 
of the planet and the lifetime of the cloud structures. 
The autocorrelation A can be computed as: 



A(L) 



J2k=o L ^ (F(kAt) F) (F((k + L)At) F) 



EL"o t{F(kAt)-F) 

(3) 

where L is the time lag in number of 3 hr points, N 
is the total number of points in the time series, and F 
is the mean flux. The autocorrelation is maximum when 
values are similar within a time lag distance and it is 
sampled according to the time resolution. This method 
was chosen before the Fourier Transform to avoid the 
harmonics of the rotational period. The duration of the 
statistically significant peaks in the autocorrelated time 
series can give us an estimation of the lifetime of the 
cloud s tructures, typicall y of around one week for Earth 
clouds (Palle et al.|2008). In the out going mid-infrared 
radiation flux, a 24 hr rotation period is clearly shown, 
a value close to the true rotation period. This rotational 
signature has two origins. First, some large regions ex- 
hibit systematically high or low brightness temperatures. 
This is the case for Indonesian and Sahara areas, as the 
former is one of the most humid and cloudy regions on 
the planet, whereas the latter is warmer and drier than 



Earth's average. Therefore, the two regions appear as 
fixed cold and hot features respectively, even on averaged 
outgoing flux maps, as we can see in Figure [ljc) and (d) 
where the TOA emission is averaged over a month and 
in Figure [He) and (f) where it is averaged over a year. A 
smaller effect comes from the diurnal cycle (the change 
of brightness temperature between day and night in a 
region), which is negligible in most locations, because of 
humidity, clouds, or ocean thermal inertia, but noticeable 
in some dry continental areas as it is shown in Figure [5j 
In fact, it is the diurnal cycle of the dry lands, which 
makes possible the detection of the rotation period for 
the case of the North-polar view during the more stable 
seasons (winter and summer), as it is shown in Section 
3.2. Although it is always the same fraction of the planet 
(the whole northern hemisphere) in the field of view, the 
change in temperature along the day in the dry continen- 
tal areas causes the rotational modulation as it is shown 
in Figure HI (violet solid line). However, an observer look- 
ing at theSouth Pole does not detect a significant rota- 
tional variability, Figure [4] (violet dash-dotted line) , not 
even for the austral summer when the effect of clouds 
would be minimized. In the Southern Hemisphere, the 
distribution of land is largely dominated by oceans whose 
high thermal inertia make diurnal temperature variabil- 
ity negligible. 

3.2. Average Rotation Light Curves 

Once the rotation period is identified, the observer can 
produce a typical rotation light curve by folding the time 
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Fig. 6. — Rotation light curves of the mid-infrared radiation emitted from the Earth. For the months of 2001 January, April, July and 
October, 0° longitude and latitudes: 90°N, 60°N, 45°N, 30°N, 0°, 30°S, 45°S, 60°S, and 90°S. 



series obtained during weeks or months over the rotation 
period. This process averages out random cloud vari- 
ability. It is clearly seen on the maps of Figures |T](c) 
(f), where the clouds disappear for longer average times 
whereas the strongest features mentioned on the previous 
section prevail. Then the observer can plot this average 
rotation light curve as a function of an arbitrary longi- 
tude (in our case the longitude that we have chosen is the 
conventional geographic longitude for commodity). The 
shape of the light curve can then reveal brighter/fainter 
areas distributed in longitude. For instance, observers 
over a latitude of 30° N would note that the brightest and 
faintest point of the light curve occur when the longitude 
0° (Sahara) or 135° (Indonesia) are respectively centered 
on his view. The results are represented in Figures [5] and 

Figure [5] represents the rotational variability with local 
hour for several latitudes (in column) and at the months 
of January (top row) and July (bottom row). Each chart 
corresponds to four observers placed at the same latitude 
referred and at the four initial longitudes O, C, WQ, EQ, 
previously defined. Each of these four planes correspond 
to a certain local hour that changes along the orbital 
period. With this information, we can compare the tem- 
perature evolution for a given planetary region along the 
day (diurnal variability). The maxima show the largest 
variation temperature along the day whereas the minima 
hardly vary. For the cases where the observers are placed 
over the poles, the planetary view does not change with 
time so the local hour of the graph is taken just as ref- 
erence of the observation time. For the North Polar case 
and for both summer and winter seasons, the minima oc- 
curs when it is UT (0 hr at 0° longitude or 12 hr at 180° 
longitude) and the maxima when it is 12 UT (0 hr at 180° 



longitude or 12 hr at 0° longitude), which illustrates the 
diurnal cycle effect previously mentioned. These hours 
coincide respectively to midnight and midday in the Sa- 
hara desert. The greatest influence of the area is noticed 
at 30° N latitude when the Sahara is in the center of the 
planetary disk. In summer or winter the diurnal vari- 
ability reaches the 2%, although the main contribution 
to the flux is that of the solar insolation along the year. 

Figure [6] represents the rotational variability for ob- 
servers placed at 0° longitude and different latitudes. 
It is important to note the temperature evolution with 
time for each geometry, implying a seasonal behavior. 
The equator shows a warm stable temperature during 
the year, whereas the summers of each hemisphere dif- 
fer being the northern summer hotter, as it is shown in 
Figure p] The maximum and minimum regions do not 
change with the seasons, except for April and October 
when the presence of the clouds can mask the signal. 

In order to check the possible source of the diurnal 
variation on the Northern Hemisphere, we have made a 
further analysis of the emitted flux by geographic lati- 
tude bands. Figure [7] shows the mid-infrared emission 
of the Earth along one orbital period. Each hemisphere 
is divided in three latitude bands of flux: 90°N-60°N 
(blue), 60°N-30°N (green), 30°N-0° (orange), 90°S-60°S 
(black), 60°S-30°S (magenta), 30°S-0° (red). We can see 
that the Northern Hemisphere exhibits a bigger seasonal 
variation than the Southern Hemisphere and this varia- 
tion is more important for mid-latitudes. As it was pre- 
viously mentioned, the greater diurnal variability comes 
from the band 30°N-0°. Besides, in this band, there is a 
decrease in flux during summer. This is due to the mi- 
gration of the Intertropical Convergence Zone to these 
latitudes, which produces large bands of high humid- 
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Fig. 7. — TOA-LW infrared emission flux of the Northern and Southern hemispheres. Total emission (top) of the Northern Hemisphere 
(black) and Southern Hemisphere (red). Latitude bands (bottom) of 90°N-60°N (blue), 60°N-30°N (green), 30°N-0° (orange), 90°S-60°S 
(black), 60°S-30°S (magenta) and 30°S-0° (red). 



ity and clouds that have lower brightness temperatures. 
As expected, the South Pole is colder than the North 
Pole. For the identification of the region that produces 
the biggest variability, we have analyzed the emission of 
six regions of the planet (Figure [8| . The chosen regions 
have an area of 30°x60° latitude~~Longitude, then fluxes 
from the same latitude band can be compared. For the 
60°N-30°N band we have chosen three regions centered 
in Europe (red), Asia (blue) and US-Canada (green). 
For the 30°N-0° band the regions are centered in the 
Sahara- Arabian deserts (red), Indonesia (blue), and the 
Caribbean- Mexico area (green) . As expected, the regions 
in the same latitude band have similar seasonal variabil- 
ities but the Sahara region emits more infrared flux and 
has greater diurnal variability. However, the emission 
received by the observer changes with the view of the 
planet. The Sahara desert lies near the limb of the plan- 
etary disk facing an observer placed over the North Pole 
and its influence on the signal is lessened by the perspec- 
tive. In Figure [9j we can see that the diurnal variability 
of the signal has a contribution of both low and mid- 
latitudes. Finally comparing Figures (5j|7] with the signal 
received, shown in Figure [3] we conclude that although 
there is a small difference Detween day and night pro- 
duced by the warm areas, the signal is dominated by the 
seasonal behavior. 

3.3. Phase Variation 

While Earth's visible flux received by a remote ob- 
server is modulated by the changing phase of the planet, 
the Earth does not present significant phases when 



the integrated infrared flux is observed: the emission 
from the nightside contri butes nearly as much as the 
emission from the dayside ( |Selsis|2004 ) . This is shown in 
Figure [3] (middle and bottom) , which represent opposite 
observers and then opposite visible phases of the planet. 
While observer (O) sees the winter midnight, the 
observer (C) sees the winter noon which shows the same 
average temperature but bigger rotational variation. 
Figure [5] shows this difference between 12 hr and hr, 
however the seasonal variation dominates. Continental 
surfaces and the boundary layer above them (roughly 
the first 1500 m of the atmosphere) experience a drop 
of temperature between the day and night. This diurnal 
cycle is insignificant above the ocean (due to the high 
thermal inertia of water) , hence over ~ 70% of the Earth 
surface. Day-night brightness temperature variations 
affect the outgoing thermal emission only in the 8-12/im 
atmospheric window and above dry continents. This 
happens either over very cold regions which in this 
case do not contribute much to the global emission, 
or over deserts which represent a small fraction of the 
Earth surface. In addition, the diurnal cycle is much 
less pronounced above the boundary layer, at altitudes 
where most of the thermal emission is emitted to space. 
This is the reason why phase-correlated variations of 
Earth brightness temperature are negligible compared 
with seasonal changes. Even for the observing latitudes 
in the 30°N-0° range, in which the Sahara diurnal 
cycle appears in the modulation, the winter midday is 
colder than the summer midnight as it is seen in Figure [5j 
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Fig. 8. — TOA-LW infrared emission flux of some continental regions in the Northern Hemisphere. In the 60°N-30°N latitude band (top): 
Europe (red), Asia (blue), and US-Canada (green). In the 30°N-0° latitude band (bottom), Sahara-Arabian (red), Indonesian (blue), and 
Caribbean-Mexico area (green). 
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Fig. 10. — Earth— Moon mid-infrared emission light curves for one planetary orbit. 45°N (top row), Equator (second row), 45°S lati- 
tudes (bottom row), in columns the signals received by the observer's placed in opposition, western quadrature, conjunction, and eastern 
quadrature at the initial time. The colors correspond to the Earth (black) and Earth— Moon system, with lowest (red) and highest (blue) 
inclination angles of the Moon's orbit according to the observer's geometry (the possible orbits are comprehended between the two). 



3.4. Earth-Moon System Light Curves 

We have modeled the mid-infrared flux of the Moon. 
Due to a very low surface thermal inertia, the tempera- 
ture map of the starlit hemisphere of the Moon can be 
calculat ed by assuming loca l radiative equilibrium at the 



surface (Lawson et al. 2007). When calculating the disk- 



integrated emission, the contribution of the dark side can 
be safely neglected due to the high temperature differ- 
ence. Thus, the flux received depends only on the phase 
of the Moon as seen by the distant observer. The am- 
plitude of the lunar phase variations depends on the el- 
evation of the observer above the lunar orbit, when we 
computed the Earth signal as a function of the latitude of 
the observer. To a given latitude of the observer, it corre- 
sponds to a range of possible elevations above the lunar 
orbit. In Figure [10] instead of calculating this eleva- 
tion consistently with the chosen observer geometry, we 
bracket the orbital light curve of the unresolved Earth- 
Moon system with the two curves obtained by adding the 
lunar signa l for the two extr eme possible elevations. As 
pointed by |Selsis| ( |2004[ ) and |Moskovitz et al | ( |2009[ ), it 
presents phase variations dominated by the Moon. We 
note that the modulation from the satellite becomes neg- 
ligible for a Moon-like satellite with 20% of the Moon 
radius, a ~ 5% of the radius of the planet. 

The two main annual variations that modulate the IR 
emission from the point-like Earth-Moon system are due 



to the seasons of the Earth and the phases of the Moon. 
These modulations present a phase shift that depends 
on the observer geometry. For some geometries, these 
two modulations are coincidental. This happens if the 
maximum of the lunar phase corresponds to Earth's an- 
nual emission maximum, which is for instance the case 
for an observer looking at northern latitudes that sees 
the Sahara at noon in July (see Figure 10 the first two 
panels of the left column). This particular observer will 
see only phase-correlated variations and may attribute 
this variability to a day-night temperature difference and 
conclude that the planet has less ocean coverage and a 
thinner atmosphere. With such particular geometry, sea- 
sonal variations could also be mistakenly attributed to 
the phases in the absence of a moon, unless the lunar 
origin of the modu lation is identified using spectroscopy 
( |Robinson||2011 b |. 



4. CONCLUSIONS 

In this paper, we have constructed a 3 hr resolution 
model of the integrated mid-infrared emission of the 
Earth over 20 years in the direction of a remote observer 
randomly located. 

The seasonal modulation dominates the variation of 
the signal. As expected, it is larger for the polar views 
because the planetary obliquity causes a bigger annual in- 
solation change for these latitudes. For equatorial views, 
the seasonal maximum occurs during the summer of the 
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Northern Hemisphere, as the latter contains large conti- 
nental masses whereas the Southern Hemisphere is dom- 
inated by the oceans. 

The rotational variability is detectable because of the 
uneven distribution of oceans and continents with geo- 
graphical longitude. The daily maximum of the mid- 
infrared flux is shown when dry large masses of land, 
such as the Sahara desert, are in the observer's field of 
view. The daily minimum appears when cloudy humid 
regions such as the Indonesian area is visible, as iced big 
zone are confined to the poles. In the polar views, the 
distribution of land does not change with time but the 
diurnal temperature variation of large continental areas 
affects the signal, allowing the detection of the rotational 
period in the North Polar case. We find that the rota- 
tional variations have an amplitude of several percent, 
which is comparable to that of the seasonal variations 
for some latitudes. It is important to remark the strong 
influence of the weather patterns, humidity and clouds 
are sometimes able to mask the 24 hr rotation period of 
the signal for several days at a time. However, this effect 
can be solved by time folding. 

It is important to point out that the Earth does not 
exhibit a significant modulation associated with phase 
variation (phase curve). This is because the integrated 
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thermal emission does not generally probe the boundary 
layer (first km of the atmosphere) where the diurnal cy- 
cle takes place. If unresolved, the Earth-Moon system 
would however present a phase variation of Lunar origin. 
A satellite of the size of the Moon would introduce a 
strong phase variability that would completely dominate 
over the planet's signal. This effect adds high complexity 
to its interpretation by photometry. 

At the light of these results it seems that future in- 
frared photometric observations of terrestrial planets can 
be useful in order to characterize their atmospheric and 
surface features. If the planet is not completely covered 
by clouds, as Venus is, the presence of strong surface 
inhomogeneities (continents) can be extracted from the 
daily variations. The seasonal cycle can also give esti- 
mates of the planet effective temperature, the variability 
in the obliquity of its orbit, and the distribution of land 
at larger scale. A further study with a Global Circulation 
Model combining Earth's emitted flux is ongoing. 
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